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Abstract 

The recent discovery of a non-zero value of the mixing angle #13 has ruled out tri- 
bimaximal mixing as the correct lepton mixing pattern generated by some discrete flavor 
symmetry (barring large next-to-leading order corrections in concrete models). In this 
work we assume that neutrinos are Majorana particles and perform a general scan of 
all finite discrete groups with order less than 1536 to obtain their predictions for lepton 
mixing angles. To our surprise, the scan of over one million groups only yields 3 interest- 
ing groups that give lepton mixing patterns which lie within 3-sigma of the current best 
global fit values. A systematic way to categorize such groups and the implications for 
flavor symmetry are discussed. 
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1 Introduction 



The origin of flavor is one of the important questions of beyond the Standard Model physics. 
All entries of the lepton mixing matrix, or better known as the Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) matrix, are of order one, with the exception of U e s. Compared to the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix whose off-diagonal entries are small, the very different 
form of the PMNS matrix seems to suggest a different origin of the two matrices. One popular 
approach to the flavor puzzle is to invoke (spontaneously broken) symmetries to describe the 
observed patterns. The leptonic mixing angles can be determined solely from flavor symmetry 
considerations (up to permutations of rows and columns of the mixing matrix). This is possible 
if the charged lepton and neutrino mass matrices exhibit the misaligned remnant symmetries 
under which charged leptons and neutrinos transform as three inequivalent singlets, as will 
be reviewed in the next chapter. 

Assuming the remnant symmetries to be part of the original symmetry group (and not a 
result of an accidental symmetry) one can then determine mixing patterns from the structure 
of discrete symmetry groups. For review on discrete flavor symmetries and their application in 
model building see [1—3]. For example the symmetry group [4-10] and S4 [11-13] can lead 
to the tri-bimaximal mixing pattern (TBM) by Harrison, Perkins and Scott [14; 15]. With 
the latest global fits results [16—18] (see Table 1) on the non-zero mixing angle #13 measured 
by DAYA BAY [19], RENO [20] and DOUBLE CHOOZ [ ], TBM is ruled out and one is 
prompted to rethink the approach to lepton flavor based on discrete groups. One possibility is 
to build models which lead to TBM on leading order and allow for large next-to-leading order 
(NLO) corrections. Here the problem usually is that quite often there are many different 
NLO corrections, which limits the predictivity of the models. Another approach is to look 
for new groups that predict a different type of leptonic mixing pattern i.e. a new starting 
point about which models could be built. In this paper, we shall follow the second route 
and therefore perform a scan of all possible finite discrete groups of the order less than 1536 
with the help of the computer algebra program GAP [22—25]. To our surprise, only three finite 
discrete groups can yield the neutrino mixing angles allowed by the experimental constraints. 

This paper is organized as follows: in Section 2 we will present the group theoretical 
procedure to obtain the PMNS matrix from a finite symmetry group. This section might be 
skipped by readers familiar with the methodology. The method of scanning through all the 
groups of order less than 1536 and the relevant results are presented in Section 3 and finally 
we conclude in Section 4. 

2 Leptonic Mixing from Remnant Symmetries 

Lepton mixing can be derived from a flavor symmetry breaking to remnant symmetries in the 
charged lepton and neutrino masses respectively. In concrete models, this is usually achieved 
via a spontaneous breaking using some scalar fields charged under this symmetry into different 
directions of flavor space. The charge assignments are chosen such that there are different 
residual symmetries in the charged lepton and neutrino sectors. The misalignment between 
the two residual symmetries generates the PMNS matrix [11—13; 26-28]. In this method, 
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AmS^ sin 2 6*i2 sin 2 # 2 3 sin 2 13 5 

[10- 5 eV 2 ] [10" 3 eV 2 ] [lO" 1 ] [KT 1 ] [lO" 2 ] [tt] 

best fit 7.62±;^ 2.55±;g| 3.20±;Jf 6.13±; 22 2.46±;|| 0.8±^ 2 
3cr range 7.12 - 8.20 2.31 - 2.74 2.7- 3.7 3.6 - 6.8 1.7- 3.3 0-2 

Table 1: Global fit of neutrino oscillation parameters (for normal ordering of neutrino masses) adapted 
from [17]. The errors of the best fit values indicate the one sigma ranges. In the global fit there are two nearly 
degenerate minima at sin 2 #23 = 0.430i",oioi see Figure 1. 



only the structure of flavor symmetry group and its remnant symmetries are assumed and we 
do not consider the breaking mechanism i.e. how the required vacuum alignment needed to 
achieve the remnant symmetries is dynamically realized. 
The PMNS matrix is defined as 

£/pMNS = V}V V (1) 

and can be determined from the unitary matrices V e and V u satisfying 

V^M e M\V* = diag(m 2 ,m 2 ,m 2 ) and V?M V V V = diag(mi, m 2 , m 3 ), (2) 

where the mass matrices are defined by C = e T M e e c + \v T M v v. We will now review how 
certain mixing patterns can be understood as a consequence of mismatched horizontal sym- 
metries acting on the charged lepton and neutrino sectors [11-13; 26-28] . Let us assume 
for this purpose that there is a (discrete) symmetry group Gf under which the left-handed 
lepton doublets L = (u, e) T transform under a faithful unitary 3-dimensional representation 
p:G f ^ GL(3,C): 

L -> p(g)L, g £ Gf. (3) 

The experimental data clearly shows (i) that all lepton masses are unequal and (ii) there is 
mixing amongst all three mass eigenstates. Therefore this symmetry cannot be a symmetry 
of the entire Lagrangian but it has to be broken to different subgroups G e and G v (with 
trivial intersection) in the charged lepton and neutrino sectors, respectively. If the fermions 
transform as 

e -> p(g e )e, v -)• p(g u )u, g e eG e ,g u eG„, (4) 

for the symmetry to hold, the mass matrices have to fulfil 

p{g e ) T M e Mtp{g e )* = M e Ml and p{g v ) T M u p{g v ) = M u . (5) 

Choosing G e or G v to be a non-abelian group would lead to a degenerate mass spectrum, 
as their representations cannot be decomposed into three inequivalent 1-dimensional repre- 
sentations of G e or G v . This scenario is not compatible with the case of three distinguished 

4 We here follow the presentation and convention in [26; 27]. 
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neutrino and charged lepton masses and we therefore restrict ourselves to the abelian case. 
We further restrict ourselves to the case of Majorana neutrinos, which implies that there 
cannot be a complex eigenvalue of the matrices p{g v ) and they therefore satisfy p(g u ) 2 = 1> 
and we can further choose det p(g u ) = 1- By further requiring three distinguishable Majorana 
neutrinos the group G v is restricted to be the Klein group Z2 x Z2. To be able to determine 
(up to permutations of rows and columns) the mixing matrix from the group structure it is 
necessary to have all neutrinos transform as inequivalent singlets of G v . The same is true for 
the charged leptons which shows that G e cannot be smaller than Z3. We can now determine 
the mixing via the unitary matrices fi e , £l u that satisfy 

^\p{9e)^e = P(g e )diag, tilp(g v )Vl v = p{g v )diag (6) 

where p(g)diag ar e diagonal unitary matrices. These conditions determine Q e , Q. v up to a 
diagonal phase matrix K eiV and permutation matrices P e , v 

It follows from Eq. (5) that up to the ambiguities of the last equation, V^„ are given by Sl e ,v 
This can be seen as 

VL T e M e M\K = ^ T eP T M e Mtp*K = pl ag ^ T e M e MtKPkia 9 

has to be diagonal (only a diagonal matrix is invariant when conjugated by an arbitrary 
phase matrix) and the phasing and permutation freedom can be used to bring it into the 
form diag(r?i,g, m^, m 2 ), and analogously for fi„. From these group theoretical considerations 
we can thus determine the PMNS matrix 

C^PMNS = (8) 

up to a permutation of rows and columns. It should not be surprising that it is not possible 
to uniquely pin down the mixing matrix, as it is not possible to predict lepton masses in this 
approach. 

Let us now try to apply this machinery to some interesting cases. We have seen that the 
smallest residual symmetry in the charged lepton sector is given by a G e = (T\T S = = Z%. 
We use a basis where the generator is given by 



p(T) = T 3 = 1 . (9) 




This matrix will be our standard 3-dimensional representation of Z3 and the notation T3 will 
be used throughout this work. It is diagonalized by 



£l\p{T)Q, e = diag(l, oj 2 , uj) and f2 e = VLt = —7= 

v3 



(l 1 1 

1 co 2 oj I . (10) 

l 1 LO UJ 2 
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and uj = e l2?r / 3 . Having fixed the basis by choosing the Zj, generator the way we just did, it is 
now essentially a question of choosing generators and studying the predicted mixing matrix. 
Let us first look at the case where there is only one generator S of G u , satisfying p(S) 2 = 1 
and det p(S) = 1: 



p{S) = 5 3 = 
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1 







-1 





V 





-1 



(11) 



Due to the degenerate eigenvalues there is a two-parameter freedom in the matrix Q. u and it 
will turn out to be useful in classifying our result later to write it as 



ni P (s)n u = dia g ( 



-1,1,-ij 



with tt v = n u U 1 z{0,5) 



(12) 
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(13) 



Obviously this does not completely fix the leptonic mixing matrix yet, as the first and third 
eigenvalues are the same and the corresponding eigenstates can be rotated into each other 
without breaking the symmetry. To completely fix the mixing matrix we have to enlarge G v 
by another generator. Let us look at the effect of adding the symmetry generator U with 




p{U) = U 3 



This fixes the value of 9 to zero, n^p(U)Clu = diag(- 
to the famous tri-bimaximal mixing (TBM) form 



(14) 



1, —1, 1), and thus the mixing matrix 



C^PMNS = ^t^-U = UuPS = 







2 J. 
3 

1 1 1 

Ve Vs V2 



(15) 



which corresponds to the mixing angles sin 2 9\2 = \ , sin 2 #23 = 5, and sin 2 #13 = 0. TBM 
pattern is predicted by the discrete group S4 = {S^^T^^U^) [11-13] with S3, T 3 and U3 
given in our example above. Until very recently, this pattern gave a good description of the 
mixing matrix and the fact that this mixing pattern can be obtained from simple symmetry 
considerations has prompted a lot of model building activity (for a general scan of models 
based on flavor groups up to the order of 100, see [29]). In light of the recent measurement 
of a non-vanishing #13 there has been interest in the physical situation where the (broken) 
flavor symmetry does not fully determine the mixing angles. For example if the residual 
symmetry in the neutrino sector G v is taken to be G u = (S) = Z2, we have seen that the 
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Figure 1: The deviations from TBM of the form (7pmns = UhpsUi3(9 , 8) generated by the angle 9 defined in 
Eq. (13). The yellow point represents TBM, the continuous lines give the deviations from TBM with the angle 
9 given by the color code in the top right corner for <5 = f f for n = 0, . . . , 5, where n = is the outermost 
parabola etc. The global fits by Forero et al.[17] up to 3-sigma are depicted areas encircled by the black lines 
respectively. One sees that the perturbations can go in the right direction, for 9 ~ .1 — .2 a satisfactory fit 
to the data can be produced. Note that the corrections to the solar angle are smaller than the corrections to 
the other angles. The black and blue points are mixing patterns that can be produced by the flavor group 
A(96), the green and red points can be derived from A(384) while the group A(1536) generates mixing angles 
represented by the pink and orange points. This will be reviewed in Section 3. 



leptonic mixing matrix is given by £/pmns = £^hps^13(0> 5) [ 6]. This mixing pattern 
leads to sin 2 6 12 > 1/3 and is sometimes known as tri-maximal pattern TM2. The result of 
a deviation from 9 = in terms of mixing angles is shown in Figure 1 where one can read 
off that a 13-rotation about an angle of 9 ~ 0.2 is required. In Section 3, we will show that 
all flavor groups with the order less than 1536 with an additional group A(6 • 16 2 ) lie on the 
parabola U PM ns = Ukp S Ui 3 (9,5 = 0). 

It should be clear from the discussion above that, a different choice of generators of the 
residual symmetry groups leads to different mixing pattern. For example, if we take the G u 
not to be the Klein group G v = (S, U\S 2 = U 2 = E;SU = US) ^Z 2 xZ 2 that leads to TBM, 
but the isomorphic group G u = (S, X\S 2 = X 2 = E; SX = XS) Z 2 x Z 2 with X = T 2 ST 
and 

/ -1 \ 

10 (16) 
V -1 ) 

with n x = nuU 13 (0 = -j,6 = ^) (17) 



p(X) = 

fixes 

tf xP (x)n x = ai a g(-i,-i,i) 
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Table 2: Groups generated by T3, 53 and U-i(n), that lead to new starting points. The series of groups A(6n 2 ) 
has been studied in [14] and the group (Zi$ x Zq) xi S3, apart from being defined by (T3, S3, f/3(9)}, is the 
group number 259 of order 648 in the SmallGroups catalogue of GAP [24]. 



which gives 



C^PMNS || = || 



X 




which corresponds to the mixing angles sin 2 #12 = \ , sin 2 #13 : 
mixing pattern is sometimes called tri-maximal mixing [37-40]. 
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3 New Starting Points 

In Section 2 we have seen that to uniquely (up to permutations of rows and columns) deter- 
mine the mixing patterns from group theoretic considerations, it is essential to have unbroken 
remnant symmetries in the charged lepton and neutrino sectors and to have enough symme- 
tries in each sector that there are three inequivalent representations in each sector. For the 
neutrino sector one needs a Klein group Z<i x Z2 for three generations of Majorana neutrinos 
as there cannot be a cyclic group of order larger than 2 in the neutrino sector as complex 
representations would forbid neutrino mass terms. Higher Zr? (n > 2) product groups are 
redundant as the Klein group is the maximal group that one could have for three generations 
of neutrinos. The condition of G v = Z<i [32-35; 41-43] can be used to determine the mixing 
angles up to a two-parameter freedom. This approach however does not allow a sharp mixing 
angles prediction but can accommodate the experimental results more easily. A group scan 
based on this assumption appeared in [28]. For the charged lepton sector we first consider a 
cyclic group with a Z 3 symmetry, and generalize it later to general finite abelian groups. 



3.1 G v — Z2 x Z21 G e — Z 3 

We first assume that the charged lepton sector has a Z 3 symmetry. The way to generate new 
mixing structures apart from TBM is now to embed these abelian symmetries in different 
ways into some larger group. Once the group and the embedding is specified, the Clebsch- 
Gordon coefficients of the specified group uniquely (up to permutation of rows and columns) 
specify the mixing pattern. 

We have performed such an analysis, considering all discrete groups of size smaller than 
1536 using the computer algebra program GAP. The total number of 1336749 finite discrete 
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Figure 2: The leptonic mixing angles (black circles) determined from our group scan up to order 1536 are 
shown. The red dots represent the mixing angles that we have determined from the generator 5*3, T3 and 
Uz{n). The red labels represent the integer n that generates the 1/3(71) matrix. The interpolating line is 
colored according to the value of 9 as defined in Eq. (13). See the main text for more detailed informations. 
We have also omitted the labeling of larger n that generates the same repeating groups or mixing angles. 
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groups (with an additional group of A(6 • 16 2 )) has been scanned, keeping only groups which 
contain the Klein group and Z3 as their subgroups. To speed up the scanning process, we 
use the Lagrange theorem to skip over groups with order that is not divisible by 4 (order of 
Klein group) and 3 (order of Z3). We then look for groups that carry at least a 3-dimensional 
irreducible representation. Furthermore we only consider 3-dimensional representations of a 
specific group if they are faithful to avoid duplicating subgroups which have been scanned. 
All the 3-dimensional representations p(U), p(S) and p(T) for generators U, S and T are 
then being recorded. Note that one should not be confused by p(S), p(U) and p(T) with 
S3, Us and T3 as 53, U3 and T3 are specific matrices defined in previous section. With 
each set of generators for the Klein group, one can determine the unitary matrix Vt u which 
simultaneously diagonalizes both p(S) and p(U) of the Klein group. The corresponding Q u will 
be subsequently multiplied with a unitary matrix f2 e that diagonalizes the p(T) matrix in the 
same representation such that a PMNS matrix can be obtained. This process is repeated for all 
the different combinations of Cl u and tl e obtained from different sets of {p(U), p(S)} and p(T) 
respectively. All the matrices with permutations of rows and columns for the corresponding 
PMNS matrices were subsequently generated and we ordered the PMNS matrices with the 
smallest 13-entry and an additional condition such that the 11-entry is larger or equal to 12- 
entry. A huge amount of duplicates is removed with this method. The flavor mixing angles 
are extracted from the remaining PMNS matrices and the results are plotted as black circles 
in Figure 2. From Figure 2 we observe that from about one million finite discrete groups that 
we have scanned, only 3 finite groups up to the order of 1535 including A(6 • 16 2 ) generate 
the mixing angles which are allowed by the experimental values up to 3a, with the strongest 
constraint coming from #13. Amazingly many of the discrete finite groups with G e = Z3 and 
G u = Z2 x Z2 that we have scanned yield a combination of leptonic mixing angles that lie 
on a parabola. However in general there exist also mixing angles that lie beyond our plotted 
region, these points will be shown in Figure 3. 

Those discrete groups that lie on the parabola can be easily presented in a simple system- 
atic way: in Section 2, we had seen that if the Z3 symmetry in the charged lepton sector is 
generated by the matrix T3 given in Eq. (9) and the Klein group is generated by the matrices 
S3 and U3 of Eqs. (11) and (14), the resulting mixing matrix is of the TBM form Eq. (15). 
Now, all new mixing patterns found in the scan can be written as 

£/pMNS = OhpsM* = \ arg(^), 5 = 0) (19) 

which is the result of the remnant symmetry T3 in the charged lepton sector and the choice 
of S3 and 



with (z) = Z n , n G N (20) 



as generators of a Klein group in the neutrino sector. One can think of ^(n) as one of the 
more generalized generator for the Klein group. Note that for any z with modulus one, we 
have [83,1/3(12)] = and ^(n) 2 = I3 and therefore the group generated by S3 and U^(n) is 
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0.3420 


0.0254 


0.3867 








0.3420 


0.0254 


0.6134 



Table 3: Mixing angles which are compatible with experimental results generated by flavor groups up to order 
1536. The group identification function in SmallGroups is not available for group with order 1536. 

always a Klein group, (83,1/3(71)) = Z 2 x Z 2 . For the group generated by T3, S3 and 1/3(71) 
to be finite, z has to be of the form given in Eq. (20), as may be seen by looking at the group 
element (Us(n)T3) 2 = diag(z, z, z* 2 ) which is of finite order n 6 N iff z n = 1. The requirement 
(z) = Z n further fixes n to be the smallest n for which z n = 1. Note that different n-th root 
of z will in general lead to different leptonic mixing angles via Eq. (19), as can be seen in 
Figure 2. The group generated by T3, S3 and 1/3(71) is always the same due to the requirement 
(z) = Z n . The names of the groups generated for n = 4, . . . , 16 can be found in Table 2 and 
the groups A(96) (n = 4) and A(384) (n = 8) have been obtained before in [27]. Note that 
not all of the groups generated in this way can be classified as A(6 • n 2 ), e.g. the group 
(Zis x Zq) x S3 . Another surprising observation is that all the groups which are restricted 
to those shown in Figure 2 give a prediction of 5cp = vr or <5cp = 0. In general other groups 
that we have scanned which lie outside the region shown in Figure 2 do give non-trivial Dirac 
CP phases. The predictions for mixing angles for all groups of size smaller or equal than 1536 
groups is presented in Figure 2. It should be clear that if one allows for groups of arbitrary 
size, the parabola depicted in Figure 2 will be densely covered. As one can see, the mixing 
patterns corresponding to n = 5, n = 9 and n = 16 give a good descriptions of the leptonic 
mixing matrix. See Table 3 for the resulting mixing angles. 

3.2 G v = Z 2 x Z 2 , \G e \ > 3 

In this section we will discuss the result of all the neutrino mixing angles scanned by relaxing 
the condition on the charged lepton subgroups. We allow G e to be any abelian groups which 
are the subgroups of the original flavor group Gj. The Klein group is kept as a subgroup of 
the remnant symmetry for the neutrino mass matrix. The scanning procedure is performed 
as in Section 3.1, with a few subtleties changed. For instance if G e consists of two or more 
generators above, the unitary matrix Q e is constructed as a matrix which diagonalizes all the 
generators of G e simultaneously. We have scanned all the discrete finite groups Gf up to 
order 511 with all the abelian subgroups G e . The result is shown in Figure 3 and we obtain 
only the neutrino mixing angles predicted by finite modular groups and their subgroups, as 
discussed in [27]. Mixing angles which are not generated by finite modular groups appear 

5 Let us mention that (Zig x Z%) xi S3 is a subgroup of A(6 • 18 2 ), however we only consider the smallest 
group that generates a given mixing pattern. 
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Figure 3: Mixing angles obtained from groups with all the abelian subgroups of G e up to order 511. For 
comparison purpose we also plotted the mixing patterns obtained from groups with G e = Z3. Only the mixing 
patterns generated by finite modular groups and their subgroups are obtained for \G e \ > 3, see Ref. [27]. The 
mixing patterns from G e = Z3 that lie outside the parabola are also generated by finite modular groups (e.g. 
A$, PSL(2, Zr)) and their subgroups. 
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however if we scan the groups with order higher than 511. For instance in Section 3.1 we have 
seen groups such as A(6 • 10 2 ) and (Z\$ x Zq) x S3 that are not generated by finite modular 
groups but give a set of mixing angles which are in good agreement with experiments. 

4 Conclusions 

We have presented an extensive scan over discrete symmetry groups that may be used to pre- 
dict leptonic mixing angles. It is assumed that the left-handed leptonic doublets are assigned 
to 3-dimensional representations and that the symmetry is broken to different mismatched 
symmetry groups G e and G v in the charged lepton and neutrino sectors, respectively. With 
the assumptions of G e = Z3 and G u = Z2 x Z2, we have scanned all groups with size less 
than 1536 with the additional group A(6 • 16 2 ) and surprisingly we found only 3 groups that 
lead to acceptable mixing patterns. In particular the groups A(6 • 10 2 ), (Zi$ x Zq) x S3 and 
A (6 • 16 2 ) were found to predict the leptonic mixing angles within 3-sigma region of global 
fits. In all the interesting cases we predicted Sep = or <5cp = 7r. These groups should be 
interesting for model building and phenomenological studies. 

We have also found a way to systematically categorize groups that generate the leptonic 
mixing angles which lie on the parabola that can be described by a 13-rotation of TBM mixing 
matrix, i.e. Eq. (19). It is possible that if our result can be extrapolated to arbitrary large 
order of discrete finite groups, with the assumption that the remnant symmetries of G e and 
G v are not broken, all of the interesting groups should generate lepton mixing patterns that 
lie on the parabola with trivial Dirac CP phase. 

In a second scan we relaxed the condition G e = Z3 and assumed G e to be any abelian 
group. In this scan up to groups of size 511 no new interesting groups were found. 

From all the discrete groups that we have scanned, only groups with relatively large order 
predict acceptable mixing angles. This might seem not economical from the model building 
perspective. However it is possible that some of the remnant symmetries are accidental, e.g. 
in the case of A± models that predict TBM, the generator U3 is not part of the group A± but 
rather represents an outer automorphism of the group A± [45]. It is possible that there exist 
some smaller groups whose automorphism group is isomorphic to the groups we have found 
(for example the smallest group whose automorphism group contains A (96) is of the order 
of 32 [ ]). This idea is interesting from group theoretical and phenomenological perspective 
and such an investigation is left for future studies. 

Acknowledgements: M.H. and K.S.L. acknowledge support by the International Max 
Planck Research School for Precision Tests of Fundamental Symmetries. 
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